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ABSTRACT 


Several  models  for  the  analysis  of  the  modulation 
of  the  ocean  surface  wave  spectrum  by  a  prescribed  surface 
current  are  reviewed.  Both  an  eigenmode  description  and 
WKB  approximation  will  be  studied  in  two  horiiontal 
dimensions.  The  results  will  be  expressed  as  a  modulation 
of  the  surface  wave  amplitude  power  spectrum. 
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1.0  Introduction 


In  this  paper  *:e  review  several  models  for  analysis 
of  the  modulation  of  the  ocean  surface  wave  spectrum  by  a 
prescribed  surface  current.  A  quasi-linear  approximation 
will  be  used.  The  free  surface  waves  will  be  treated  in 
the  linear  approximation  and  the  coupling  to  the  surface 
current  will  be  proportional  to  the  product  of  the  surface 
wave  and  surface  current  amplitudes.  Both  a  modal  descrip¬ 
tion  and  the  WKB  approximation  will  be  studied.  The  results 
will  be  expressed  as  a  modulation  of  the  surface  wave  ampli¬ 
tude  power  spectrum. 

The  notation  of  Parts  1^  and  II  ^  will  generally 


be  followed. 


2.0  Description  of  the  Problem 


The  model  system  we  consider  in  this  paper  is  that  of 
an  irrotational ,  incompressible  ocean.  We  treat  the  surface 
in  two  horizontal  dimensions  with  the  undisturbed  ocean  sur¬ 
face  chosen  to  be  the  z=0  plane  of  a  rectangular  coordinate 
system  with  z-axis  directed  upward.  The  depth  of  the  water 
is  assumed  to  be  very  large  compared  with  all  wavelengths 
of  interest.  The  assumption  of  irrotational  ty  allows  us  to 
use  a  velocity  potential  description  of  the  surface  waves 
[<j>(r,t)j  and  the  vertical  displacement  of  the  water  surface 
from  equilibrium  is  written  as  :;{r,t).  The  prescribed  surface 
current  is  assumed  to  be  parallel  to  the  x  axis  and  a  super¬ 
position  of  modes  of  the  form 


U(r,t) 


U  cos(K^) 
K 


where 

£  ®  x  -  Cjt 


(2.1) 


(2.2) 


with  Cj  a  phase  velocity,  K  representing  a  set  of  wave  numbers 

and  the  U  a  set  of  amplitudes. 

K 

The  linearized  equations  to  determine  4>  and  £  were 

(3) 

obtained  for  one-dimensional  waves  by  Zachariasen  and 
Milder ^  and  also  in  Part  II.  These  can  be  generalized  to 


two  dimensions;  the  resulting  equations  are 


_  +  £.vs  $  +  g  ?  =  0, 

If  +  U-Vs  ?  +  5  va=®8*  •  ,2-3 

Here  g  is  the  acceleration  of  gravity,  <j>  is  the  velocity 
potential  at  the  surface,  and  z,  is  the  vertical  displacment 
of  the  surface  from  the  plane  z=0.  The  operator  Vg  is  the 
gradient  in  the  (x,y)  plane  and 


®s 


Following  the  notation  of  Part  I,  we  write 


<f>  (r 


'*=>  -  Z  ^  ekZ  b,+»  +  b,+>* 


/2  k  \  k  -k 


and  for  the  surface  displacement  from  equilibrium 


S(r,t)  = 


I 


.  b<+>* 

/2  k  jt  -k 


where  v.  =  /g/k ,  the  b^+^  are  expansion  amplitudes,  Qnd  the 
K  k 

sum  runs  over  s  irf ace  wave  number  vectors  k.  Substitution  of 


(2.4) 


(2.5) 


Eqs.  (2.4)  and  (2.5)  into  (2.3)  leads  to  the  equation  for  the 


i  £i+>  -  „k  bi+' 

k  K  k 


kuK 

k  -K  -ifi„t 
x  „  K 

kx  +  "k  1 

4 

|S-S| 

kx'K  “|£-K|J 

k  +K  in„t 

+  -X -  e  K 

Iic+K| 


k  ' 
x 


(a). 


_  T  _  b(  +  ) 

Lkx+K  "  “\t+K\\  k+K 


(2.6) 


where  we  have  seen  in  Part  II  that  the  coupling  between 

and  b^*  ( =b  ^ ~ M  is  very  weak  and  therefore  ignored  and 
it  -t  '  it  ' 


=  Kc  and  id.  =  "V  g  |  Jc  | 
K 


(2.7) 


Continuing  to  follow  the  notation  of  Part  I,  we  set 
the  eigenmode  expansion  amplitudes  to  be  of  the  form 


b<+>  =  e'l“kt  q  a-15 

£  £ 


and  define  the  complex  amplitude 

i  (&*r-<Dkt) 


(2.8) 


Z (r,t) 


'■*  “  wk~'  ,  -l 

q£ 


(2.9) 


In  terms  of  this  complex  amplitude  we  may  write  the  surface 
displacement  as. 


C (r,t)  =  -  Im 


jz  (r,t)} 


=  -  I  Z  (r,  t)  -  Z 


(r,t)) 


(2.10) 


4 


Until  now  our  discussion  has  been  concerned  with  the 
mechanical  interaction  of  waves  and  the  resulting  surface 
displacement.  Suppose  we  shift  our  focus  and  determine  the 
spectral  development  induced  by  the  interaction  of  the  sur¬ 
face  current  and  surface  gravity  waves.  Let  us  suppose  that 
at  a  given  time  t  we  observe  q  to  have  the  value  Q  .  A 

O  V  r*- 

k  k 

series  of  such  observations  at  time  intervals  long  compared 
to  correlation  times  will  lead  to  an  ensemble  of  values  for 
the  q  's.  We  shall  assume  the  different  Q  's  to  be  uncorrelated 

H  Jc 

so  that  for  an  ensemble  average  denoted  by  the  brackets  <•••> 
we  have  the  relations, 

<  Q  >  =  <  Q  Q  >  =  0  » 

£•  k  t 


<Q,  Q*>  =  <|Qj2>  6  .  (2.11) 

&  X  k  kl 

We  can  use  these  considerations  to  construct  the  cor¬ 
relation  function  for  the  surface  displacement  C  between  two 
points  separated  by  a  distance  x  as  follows, 


C(x)  =  <c (r , t)  c(r+x,t)> 


=  ^\[p(&)  +  P(-£)}  ex^': 


(2.12) 
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where  P(ic)  is  the  power  spectrum  of  the  complex  amplitude  Z. 
Using  the  expression  for  the  surface  displacement  in  terms 
of  the  complex  amplitude,  i.e.,  Eq.  (2.10),  in  Eq.  (2.12) 
yields  the  expression, 


P(k)  = 


1  ,2 

- d  x 

2  (2tt )  J 


e  ik’x  <Z  (r,t)  Z  (r+x , t) > 


=  <|Qj2>/2k2  , 
k 


(2.13) 


which  is  the  power  spectrum  in  terms  of  the  measured  slope 
variables . 

The  corresponding  power  spectrum  of  the  surface  dis¬ 
placement  £  is  defined  by  Phillips ^  as 


H'(k) 


=  >5 


P  (k) 


+  P( 


(2.14) 


The  normalization  of  the  functions  T ( k )  and  P(k)  are  such  that 


P  (k)  =  <;2> 
k 


(2.15) 
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3 . 0  Analysis  of  the  Coupling 

In  the  remainder  of  this  paper  will  assume  that  there 
is  a  single  nondispersive  internal  wave  which  gives  rise  to 
the  surface  current.  We  write  Eq.  (2.1)  as 


U  -  i  UQ  cosK£  (3.1) 

where  the  wavelength  of  the  internal  wave  ( 2tt/K )  ,  i.e.,  sur¬ 
face  current,  is  very  much  larger  than  the  surface  waves  of 
interest.  Thus,  we  have  the  relation 

k  >>  K  ,  all  k  of  interest.  (3.2) 


To  simplify  our  two-dimensional  problem  we  parameterize 
the  dependence  of  the  surface  wavenumber  in  the  direction 
orthogonal  to  the  current  by  introducing 

p  E  3*k  (3.3) 

into  Eq.  (2.6).  With  p  as  a  fi<ed  parameter  in  Eq.  (2.6), 
the  component  of  k  parallel  to  th>  current  U  is  written  as. 


7 


Eq.  (3.4)  is  reasonable  because  of  condition  (3.1);  i.e., 
we  can  establish  the  periodic  boundary  conditions  for  the 
surface  gravity  waves  over  an  arbitrary  patch  of  ocean. 

Implementing  Eq.  (3.4)  allows  us  to  replace  the  eigen¬ 
mode  amplitudes  b^  by  a  discrete  set  of  quantities  B(n): 


-inn  c 

e 


B  (n) 


(3.5) 


where 

n  =  KcT 


(3.6) 


is  the  frequency  of  an  internal  wave  of  wavenumber  k 
and  phase  velocity  Cj.  Substituting  the  discrete  variables 
defined  by  Eq.  (3.5)  into  the  interaction  equation  ^Eq.  (2.6) 
results  in 


i  B  (n)  -  En  B  (n)  =  f(n)  B(n+1)  +  g(n)  B(n-l)  (3.7) 


vhere 


KU 


g(n)  -  n  -j 


o  n-1 


V(n-l)  2  +  (p  /  K)  2 


K  ♦  Js_ ) 

l"-1  “n-l  /  ' 


(3.8) 


KU 


f(n)  =  n 


n+1 


V(n+1)  2  + 


(n+1 )  + (p/K) 


!.) 


n 


n+1  u) 


n+1 


(3.9) 
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and 


E  =  u>_  -  nft 

n  n 


with  u>  =  a)  .  A  further  transformation  on  Eq.  (3.' 

n  IS 

it  into  the  form, 


i  |i  =  w 


where  H  is  Hermitian.  Indeed,  the  transformation 


B  (n)  =  yn  ip  (n) 


leads  to  the  equation,  using  Eq.  (3.7), 


1  *<»>  -  En  *<">  -  Vn , n+1  *(n+1)  +  Vn,n-1  *  " 


where  the  matrix  elements  are  given  by 


and 


Vn,n+1=  f(n)  Yn+1/Yn 


Vn,n-1  =  5(n)  Yn-1/Yn  ' 


The  matrix  V  will  be  symmetric  if  we  choose 


(3.10) 

)  puts 

(3.11) 

-1)  (3.1'’ 

(3.13) 


9 


Yn~]/Tn  =  f  (n-1)  /g  (n) )  ** 


(3.1 


so  that 

Vn+1  =  /£(n)  9<n+1> 

and  (3.1 

Vn,n-1  =  9(n)  ' 

Eq.  (3.12)  now  has  the  foim  of  a  "Schrodinger  equation", 


dip 

3t 


(K  +  V)  ip 


(3.1 


where 


K  , 
nn ' 


=  E 


n  nn 


Eigensolutions  to  Eq.  (3.16)  are 

ip  ,  »  -iE^t  ,  ,  . 

r  (n)  =  e  (n )  , 


(3. 


where  X  labels  the  eigenfunctions  ^  and  eigenfrequencies 
EX .  Substitution  of  (3.17)  into  (3.16)  leads  to  the 
equation 

<EA-En)  *x(n>  =  Vnn+1  *x(n+l)  +  Vn<n.x  *x(n-l)  . 

Because  of  physical  limitations  on  our  model,  as  well  as 
mathematical  limitations  on  our  equations,  it  is  useful  to 


T 


truncate  the  set  of  Eq.  (3.18)  at  some  maximum  value  of  n, 
sav  n  .  We  suppose  that  n  „  is  chosen  sufficiently  large 
that  for  time  intervals  of  interest  this  truncation  does  not 
affect  our  conclusions.  Convenient  boundary  conditions  are 
to  be  imposed  on  the  ,  such  as 


<K  (n  .  )  =  iK  (n  )  =  0  , 
min'  y X  max 


(3.19) 


where  n  .  is  an  assumed  minimum  value  of  n. 
min 

The  eigenfunctions  ii>^(n)  are  supposed  to  be  so  chosen 
that  they  satisfy  the  ortho-normality  relations 


(n)  0>x’(n)  =  <5X  f  x  i  ' 
Vn>  =  Sn,n'  • 


(3.20) 


To  discuss  the  solutions  of  Eq.  (3.16)  several  para¬ 
meter  regimes  must  be  recognized.  There  is,  first,  the 
condition  of  "resonance",  defined  by  the  relation 


where  c^  =  (g/k)  /2  is  the  group  velocity  of  the  surface  waves. 
The  condition  (3.21)  determines  a  value  of  k  ,  say  k  =  NK, 
at  which  resonance  occurs.  Condition  (3.2)  implies  that 


11 


N  >>  1,  so  we  may  conveniently  take  N  to  be  an  integer. 
Eq.  (3.21)  can  be  rewritten  in  the  more  convenient  form 


cT  =  c  COS0 

I  g 


where 


COS0 


+  N2K2  ] 


( 


( 


tha1-  is,  0  is  the  angle  between  the  surface  current  and 
For  values  of  the  svu  :ace  wa/enumbers  near  resonance,  i. 
n  as  N ,  we  may  express  the  diagonal  matrix  element  in  Eq . 
by  the  expansion  about  n  =  N 


En  * 


+  (n-N) 


3E 
_ n 

3n 


n=N 


(n-N) 
2  ! 


32E 


n 


3n‘ 


n=N 


Using  Eq.  (3.22),  we  may  re-write  this  as 


Er  to  en  -  a  (n-N)  2  +  .  .  . 


where  a  is 


(3 


a 


Kc 


4N 


M 


3  cos  9 


( 


.22) 


.  23) 


(3.16) 


.  24a) 


.  24b) 
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From  the  form  of  a  above,  we  anticipate  that  anomalous 


resonance  effects  will  occur  for 

2  2 
cos  0  *  | 

or 

e  -  35°  .  (3.2: 

We  shall  see  that  for  6  <  35°  modes  with  n  near  N  tend  to 
be  "trapped"  in  this  neighborhood,  whereas  when  0  >  35°, 
modes  originally  near  N  tend  to  spread  indefinitely  away. 

When  we  are  interested  in  mode  numbers  n  near  n  =  M, 
we  can  use  Eqs.(3.14)  and  (3.15)  to  write 

Vntl  -  Vn±l,n  “  (Y)  /HTHTTT  ['  +<^)]  '  ,3'2' 

and 

Yn  =  (M/n)*5  exp  £(n-M)  ^  cos20  +  l^/ (2M)j  ^1  •e®]  '  ” 1 

where  we  have  specified  that  ym  =  1. 

When  the  coupling  V  is  sufficiently  weak,  we  have  the 


first  order  perturbation  solution  to  Eq.  (3.18): 


1 


We  anticipate  that  this  will  be  a  useful  approximate 
solution  when 


X  ,  X  + 1 
EX+1~EX 


<<  1 


(3.29) 


This  condition  is  most  stringent  at  X  =  N,  where  we  intro¬ 
duce  the  parameter 


V 

N ,  N+l 
EN+1_EN 


2U  N 


(3 


2n 
cos  0 


-2|] 


(3.30) 


When  S  and  (all  X)  are  sufficiently  small,  we  can 
use  perturbation  theory.  When  is  not  small  and  when  X  is 
not  near  a  resonant  mode,  the  WKB  method  is  probably  the 
simplest  to  discuss  the  surface  wave  modulation.  The  WKB 
method  fails  near  resonance,  so  for  S ,  ,  large  special  tech¬ 
niques  are  required. 

For  convenient  reference,  in  Fig.  (1)  we  show  lEn+^_Enl 
as  a  function  of  (n-N)  for  the  case  that  0=0  and  that 


Cj  =  1  m/sec 

-2  -1 
K  =  10  m  -1 

N  =  250  cos0 
2 

k  =  2 . 5  cos  9 


14 


sec 


The  quantity  (E^-En)  is  also  shown  in  Fig.  (1)  for  these 
same  parameters. 

To  discuss  the  resonant  strong  coupling  (S  >>  1) 
solution,  we  consider  mode  numbers  near  n  •=  M  »=  N  and  define 
an  integer  vQ  as  most  closely  approximating  the  relation 


v  =  /S  . 
o 


(3.31) 


For  n-M  v  we  set  E„  =  E„  and 
o  n  M 


V  V  =  U  KM/2 

n  ,  n±  1  _  o  o 


(3.22) 


in  Eq.  (3.18),  which  becomes  simply 


<PX  (n)  =  H  [  4>x  (n+1)  +  (n-1 )]  , 


where 


(3.33) 


ex  =  (EA-EM)/(2Vo)  . 


(3.34) 


The  boundary  conditcns  (3.19)  are  now  taken  as 


*x (M  ±  VQ)  =  0  . 


(3.35) 


This  choice  of  boundary  condition  wilJ  not  limit  the  validity 
of  our  results  for  initial  modes  having  n-values  well  within 


the  interval  M-v  <  n  <  M+v  ard  times  short  enough  that 

o  o 

modss  are  not  excited  near  n  =  (N±v  )  or  (M±v  ). 

o  o 

The  eigenfunctions  of  Eqs.  (3. 33) -  (? . 35)  are  then 


(n) 


X  (n 


-  M  +  v 

o 


e 


X 


cos 


n  «  M  -  v  +  1,  .  . . ,  M  +  -  1  , 


X  =  1,  2,  ...  2vQ-l  .  (3.36) 

Equation  (3.33)  can  be  generalized  if  we  use  (3.16) 

in  (3.18),  but  continue  to  set  =  E„.  The  resulting 

n  m 

equation  is 

(EX-EH)  (n)  = 


If  we  define 


4>^(s)  = 


■  I  *- 


/n  ^ (n)  , 


n 


there  results  the  equation 


vx  »x  -  E  [(1+s2>  h]  ' 


(3.37) 


(3.38) 
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We  now  have  VQ  =  UQ  KN/2,  since  we  have  taker  M  =  N. 
deriving  Eq.  (3.41)  we  have  neglected  a  term  of  order 
compared  with  (n)  and  have  neglected  derivatives  of  higher 
order  than  the  second. 

Equation  (3.41)  is  the  "harmonic  oscillator"  equation 
with  "spring  constant"  a.  For  cos9  >  Jjj  ^Eq.  (3.24)}  a 
is  positive  and  the  "restoring  force"  limits  the  spreading  of 
modes.  For  cos0  <  jj]^  the  force  is  "anti-restoring"  and 
indefinite  spreading  in  mode-space  occurs. 


In 
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If  we  define  the  quantities, 


v  =  (n-N)/S 


h=-[2- 


(EX  -  EN) (Vc  l'1  /SH 


(3.42) 


and  (v)  =  ^  (n) 


\-hen  Eq.  (3.41)  can  be  written  as 


2 

d^F 


dv 


I  +  -  v2)fx  *  0  • 


(3.43) 


This  has  the  solutions, 


=  2X  +  1  ,  X  =  0,  1,  2,  ... 


Y^(v)  =  S**  2*  X!  /if  (v)  exp[-v2/2] 


(3.44) 


(n) 


where  H.  is  the  Hermite  polynomial 


H,  (v)  =  0  . 
X  2tti 


2vz-z  dz 

e  T+T 

z 


(3.45) 
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and  the  contour  of  integration  is  a  circle  about  z  =  0 
The  normalization  of  Eq.  (3.44)  is  so  chosen  that 


dn  \1>X  in)  ,  (n)  =  6^ 


and  the  eigenvalues  of  E  are  of  the  form 


EX  =  EN  +  VQ  (2  -  S_>5)  -  a6o) 


(3.46) 


vrhere 


=  2Vq  S 


-k  fcT  U 

*5  _  I  o 


[  2 


(3  cos  9-2) 


]'■  • 


(3.4?) 


(7  1 

The  above  result  is  similar  to  that  obtained  by  Rosenbluth' 
for  the  one  dimensional  impulsively  applied  internal  wave; 
the  eigenvalues  being  identical  for  9=0. 


mumm 


so  (4.4)  becomes  now  cosQ  =  MK/k 


M 


i  (kMx-wMt) 

z  =  g  M  — —  q(m)  g(m) 


M 


(M)  _ 


• 

1  -  i  - 

cos  ^0  9  1 

> 

V  ei(n-M)Kt  t(n.t) 

4MK  9xJ 

.4-  J 

(4.5) 


If  we  expand  in  terms  of  the  eigenfunctions  of  Fq.  (3.18) 
and  use  (4.2) ,  we  obtain 


ip  (n,  t)  = 

/2 


(M) 


Z 


e“lE  ^(n)  4>x  (M) 


M  u  X 


(4.6) 


Evaluation  of  (4.5)  when  we  can  use  the  perturbation 
solution  (3.28)  is  straightforward.  There  results 


,  4  U  KM 

g(m)  a  !  .  o 


1  1  2,.,1  iKC 

,1  +  2S  11  +  2  cos  01  e 


1  -  e 


■i'Vr'n11! 


-  -  as  a  +  i^20))e'lK’ 

x  -  e‘i(EM-rEM)t‘J/(EM.1-V  ■ 

On  expanding  Eq.  (4.7)  in  t  near  t  =  0 ,  we  obtain 


(4.7) 
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1 1(1  +  5 


cos  0)/M  sin (K£) 


UoKMt2) 

-5—)  ^cos(K5) 


For  9-0  and  M  =  N,  this  agrees  with  the  results  of 
Zachariasen  when  an  error  in  his  work  is  corrected. 

To  obtain  a  more  accurate  description  near  resonance 
we  write  (K  +  V)  in  Eq.  (3.16)  in  the  form 


K  +  v  =  Ho  +  n  , 


V  =  Em  *(n)  +  Vo  ['Mn+1)  +  I|/(n-l)]  , 

h1*’  *  (W  *(n)  +  Vn,n+1  *(n+1)  +  Vn , n-1  ' 


where 


Vn,n±l=Vo 


We  next  take 


iMn)  =  4>°  (n)  +  i|^(n)  , 


treat  and  ^  as  small,  and  specify  that 


(4. 


In  anticipation  of  the  evaluation  of  (4.5),  we  use  (4.15) 
and  find  that 


lvK^  4>°  (n)  =  exp|-iU0KMt  cos  (KOJ  •  (4.17) 

v 

If  we  neglect  the  term  of  0(1/N)  and  substitute  the 
above  into  (4.5)  we  see  that  the  effect  of  the  surface 
current  is  to  modulate  the  phase  velocity  of  the  surface 
wave.  This  velocity  is 


+  U. 


cos0  cos(K£) 


(4.13) 


a  result  that  could  have  been  deduced  from  elementary 
considerations . 

To  obtain  corrections  ofO(l/N)  we  must  integrate 
Eq.  (4.13).  On  setting 


1 

^(n)  =  e 


0(M)  1 
-9 -  $i(n) 


^  YM 


(4.19) 


we  obtain 


$ 


1 


(M±l) 


S{ 


=  <?it/(2N)  + 


<EM±rV 


(4.20) 


for  small  t.  Evaluation  of  (4.5)  finally  gives  us 
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1  + 


(M)  ~ 


V  t 


—  ^  cos20  sin(KC) 
N  l 


T  -i2V  t  cos(KU 
o 
e 


f  2  32E 

+  V  I  —  sin  (K^ )  —  -p~  2 

°|^N  2  aM2 


cos  (Kf, ) 


(4.21) 


This  agrees  with  the  Zachariasen  expression  (4.8)  for  short 
times.  Since  the  second  term  above  has  been  evaluated  only 
for  small  t,  we  cannot  use  (4.21)  for  late  times.  From 
Eq.  (4.15)  we  see  that  the  probability  that  mode  number 
n  is  excited  is  proportional  to  [jy  ^2v0t^]  •  Thus,  when 


v 


the  neglect  of  (En~EM)  in  Eq.  (3.33)  is  no  longer  valid  and 
we  must  use  the  Rosenbluth  equation  (3.41)  to  study  the 
interaction . 

Before  doing  this,  it  is  instructive  to  construct 
the  wave  amplitude  Z  using  the  function  4>x  of  Eq.  (3.40). 
After  some  simplification,  and  with  the  identification  tnat 


s  =  e 


iK  F, 


we  obtain 


'lEMt 

Z  =  e  (/2  YM/kM 


)  ei(py+Kr)  [\  +  |  (  cos20-2)/(4M)] 


( s  +  i  “■  m)|  r  (s,t) , 


(4.22) 
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where 


T(s 


"  ■[? 
■/ 


e-lE^t  +x (s) 


iEMt 


dv  C  (v)  v(e-ct) 
5  e 

1+s 


(4, 


Here  we  have  replaced  the  sum  over  X  by  an  integral  over 


v  =  2(ea-em)/(uo:o 


(4. 


and  ha^e  set 


0  =  tan-1 (s)  ,  a  =  itUQK/2  .  (4. 

The  coefficient  =  C(v)  is  chosen  to  satisfy  the  boundary 
condition  that 

r  (s,o)  =  sM_1  Q(m)/(/lYM)  .  (4. 

Using  (4.26)  wa  obtain 

T(s,t)  =  tanM_1  ( 8  —a )  tan2(0-a)J  Q(M)  (4> 

(1  +  tan20)  /2  ym 

On  evaluating  (4.22)  we  obtain  the  envelope  function 


(l-o2)  £l  +  o|^  cos20-l]  sin(K£)/ 

1 1  +  o2  -  2a  sin ( )  ik  +  o2  -  2o  sin(K^)) 


(4.28) 


Here 


o  =  tanh  (UQKt/2)  . 


(4.29) 


For  small  t,  o  «  VQt/N  and  Eq.  (4.28)  is  in  agreement  with 

the  linear  t-dependent  terms  of  Eq.  (4.21).  Eq.  (4.28)  is 

singular  for  a  =  1,  sinK£  =  1.  Including  the  effect  of  the 

(E  -E  )  term  in  (3.18)  would  remove  this  singularity 
n  M 

For  a  <<  1  we  can  simplify  Eq.  (4.27)  to  obtain 


V  t  , 

1  +  (1+^cos  9)  sin(KC) 

N 


exp 


£-2iVQt  cos  (K£ )j  , 


(4 


which  modifies  Eq .  (4.21)  with  a  plausible  phase  factor 

2  2 

correction  and,  of  course,  omits  the  term  involving  9  EM/9M  . 
Returning  to  the  Rosenbluth  equation  (3.41),  the 
general  form  of  >Mn,t)  is 


<Mn,  t) 


cx  n(n) 

x 


(4.30) 


where  we  have  dropped  a  phase  factor, 


and  5u)  and  ip ^  ( n )  are  given  by  Eqs.  (3.42)  to  (3.47).  The 


29) 
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(4.31) 


A  special  case  of  some  interest  is  that  of  the  "non¬ 
spreading"  wave  packet,  with 

A  =  1,  A'  =  1  , 

'  «  v  e2i6a)t  (4.35) 


In  this  case 


1 1|»  (n,  t)  |  2 


(M-N) 


cos  (6wt ) 


(4.36) 


The  most  probable  value  of  n  is 

n  =  N  +  (M-N)  cos(5ojt)  (4.37) 

Returning  to  Eq.  (4.33),  we  re-write  this  as 
ijj(n,o)~  exp  (n-M)  2/ (2AsSj 

To  satisfy  our  boundary  condition  (4.2)  we  evidently  require 


that 


2AS*5  <<  1  .  (4.38) 

Reference  to  Eq.  (4.34)  shows  then  that  when 

2i6wt  . 
e  =  -1 

| A' |  is  very  large  and  we  encounter  Rosenbluth's  pile-up 
near  cosK£  =  1.  This  "pile-up"  is  greatest  at 
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(4.39) 


t  =  T 


p  =  2  6  oj 


Rosenbluth's  estimate  gives 


(M)  ^ 


G S  1  +  R(M),  r^>  >  o 


(M) 


=  1, 


R(M)  <  0 


(4.40) 


{M)  =  sin[2.8  S1/8  sin(K£)1 


sin (K£) 


-1 


(4.41) 


On  setting  M  =  N  in  Eq.  (4.9),  we  model  G(M)  as  follows; 


g(M)  ~ 


1  + 


1 _ 2, 


(1  +  ±  cos  0) 

T  i  N  -  sin  (K£) 


tr*  ’  2 


cos(KC)  /  _t\  2  .  t  \ 2 

s  N  ^ 


.  ty  (m) 

R  )  exP  I  "iMKU(C)t| 


J-iMKU(C)tj  , 


tf  T  (4.42) 
P 


where 


T V  =  (UoKN/2) 


-1 


(4.43) 


7T 

and  Tp  =  4  xv  s  ,  an  alternate  form  for  Eq.  (4.39). 

For  applications  it  will  be  desirable  to  re-label 
the  envelope  function  (4.5)  with  the  wave  number  k  =  (MK,p) , 
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setting 


(4.44) 


For  times  much  larger  than  Tp  ,  phase  mixing 
substantially  reduces  the  resonant  modulation  of  the 
envelope  function.  We  note  that  for  t  =  x  ,  the  variable 

r 

(4.29)  is 

a  =  tanh  |  ^  S^/N |  , 

which  we  anticipate  to  be  small  compared  with  unity.  This 
seems  to  justify  the  linear  approximation  in  (4.42'. 
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5.0  The  WKB  Approximation 

For  a  mode  number  M  sufficiently  far  from  resonance 
that  we  can  set 


EM+1 


Sem 
p  —  M 

EM“ 


Tm" 


9 


the  parameter  (3.29)  is 


SM  * 


U  M 
o 


2 (cl"cg  cos0) 


(5.1) 


If  this  is  greater  than  unity  the  perturbation  solution  (3.28) 
cannot  be  used,  but  the  WKB  approximation  can  be.  To  develop 
this  we  re-write  Eq.  (2.3)  in  terms  of  t,  £  -  x  -  Cjt  ,  and  y: 


H  +  V  If  *  *  -  0  . 


n 

at 


+  v  If  +  !  at=®s  * 


dU 


(5.2) 


Here 


V  =  U  -  Cj  .  (5.3) 

Next,  we  write 
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A 

Y 


A(rj  ei(py-ut) 


5 

ij  q(5'  )<^' 


C 


a(0 


ei  (py-wt) 


i/  q  a  ' )  dc 1 

e 


(5.4) 


As  usual,  we  consider  q,  A  and  a  to  be  slowly  varying 
functions  and  neglect  higher  than  first  derivatives  of  these 
as  well  as  products  of  derivatives. 

/  Q  \ 

Treating  A'  and  q'  as  small,  we  then  work  out  the  relation 

®s*S 
where 

2  2  \ 

k  =  (p  +  q  )  • 


kA  -  i  £  A' 


-  i  I1  -  l£]  V1 


i  (py-to  t) 


if  qd C 


(5.5) 


Substitution  of  Eqs.  (5.5)  and  (5.4)  into  Eq.  (5.2) 

is  straightforward  and  A  can  be  eliminated  from  the  coupled 

T  da  ^k 

equations.  There  results  here  a'  =  g|-  ,  etc.,  and  cg  = 


r 

j(u-qV)2  -  u^la  +  ik'  ^4  cg  (V  +  cg  cos9)|  a' 


+  a  j  2  cg  V '  +  2  c^  V  +  4  cg  c^  cosG  +  2  cg  sin20^ 


=  0 


(5.6) 
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The  first  term,  involving  no  derivatives,  gives  the  usual 
dispersion  relation 

a)  -  q  V  =  u)j^  =  */gk  (5.7) 

for  waves  traveling  in  the  positive  x-direction.  The  terms 
involving  derivatives  may  be  integrated  to  yield 


"\Jc  g(V  +  c  cose)  a(£)  =  const. 


(5.8) 


To  first  order  in  U,  Eq.  (5.8)  becomes 


a(U  S  a 


1  +  x  cose  [c.  -2c  cose  +  2  c  sin^6|/(c  -c  cose] 
4  Li  g  g  Jig 


=  a  G(Mi  , 
o 


(5.9) 


where  again  M  labels  the  mode. 

For  sufficiently  short  wavelengths,  we  may  neglect 

c  compared  to  c  above.  Then 
y 

a(5)  =  aQ[  1  +  Vcj8-]  '  (5'101 

The  short  wavelength  will  be  subject  to  modulation  effect 
due  to  interaction  with  long  wavelength  waves.  Let  l  be 
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a  wave  number  for  which  2  k. 


In  the  linear  approximation 


we  have 


"  V2  q2  cos 


=  Qi  G ^  cosU-x-^t) 


(5.  1 1  ) 


Here  V  =  (.i  .  and  G„  is  the  envelope  modulation  function 

IJCi.  c.  W  ^  $  /  )l.  * 

for  the  wave  t.  Then,  in  Eq.  (5.8)  we  have 


V  =  U,  (Grl)  -  (VrUv) 


(5.12) 


where 


U  =  V  Qn  cos  (8.  •  x-m,,  t ) 
u  V  2  v  i  - -  8 


(5.13) 


Linearizing  the  relation  (5.8)  in  (G^-l),  and  again  assuming 

dm.  ~ 

<<  V„  ,  we  find 


'g  dk 


A  o  cos  (  2  •  x-m  „  t ) 

k •  2,  w 2  - 2  u  cosh 

,  i  + 


a(f  )  =  aQ  \  1  +  4  1  -  Qg  cos(2;x-M(t)  '"X 


(5.14) 

Here  we  have  added  the  direct  interaction  given  by  Eq .  (5.10). 

On  squaring  and  performing  an  ensemble  average,  we 

obtain 
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<|a2(£)|>=  <|ao|2>|l  +  V  Ql  Re(Gc1]  +  -r*1  c~  '  • 

L  ~  1 

Here  we  have  neglected  higher  moments  than  the  second  of  the 
slope  function  .  If  we  sum  over  £  and  introduce  the 
spectral  function  (2.12),  this  becomes 


<|a2(?)  |>  5  <|a2|>  I  Gk  |  2 


-  <u0r>  (i  * 


I 


o  ~  ~  cosfi  U 

V  (~  ~)  P{~)  Re{Gl’1)  +  c“- 


(5.15) 
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